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Abstract 

In this paper we study second order stationary Mean Field Game systems under density 
constraints on a bounded domain O c We show the existence of weak solutions for 
power-like Hamiltonians with arbitrary order of growth. Our strategy is a variational 
one, i.e. we obtain the Mean Field Game system as the optimality condition of a convex 
optimization problem, which has a solution. When the Hamiltonian has a growth of 
order q' e]l, d/((i - 1)[, the solution of the optimization problem is continuous which 
implies that the problem constraints are qualified. Using this fact and the computation 
of the sub differential of a convex functional introduced by Benamou-Brenier (see [1]), we 
prove the existence of a solution of the MFG system. In the case where the Hamiltonian 
has a growth of order q' > dl(d- 1), the previous arguments do not apply and we prove 
the existence by means of an approximation argument. 


Resume 

Dans ce papier on etudie des systemes de jeux a champ moyen sous contrainte de densite 
sur un domaine borne D c R'^. On demontre I’existence de solutions faibles pour des 
hamiltoniens de type puissance avec ordre de croissance arbitraire. Notre strategie est 
variationnelle, on obtient le systeme de jeux a champ moyen comme condition d’optimalite 
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d’un probleme convexe, lequel a une solution. Quand I’hamiltonien a un ordre de crois- 
sance q' € ]1, d/la solution du probleme d’optimisation est continue, ce qui implique 
que les contraintes du probleme sont qualifiees. En utilisant cette propriete et le calcul du 
sous-differentiel d’une fonctionnelle convexe introduite par Benamou-Brenier (voir [1]), 
on demontre I’existence d’une solution du systeme MEG. Dans les cas on I’hamiltonien a 
un ordre de croissance q' > dj^d- 1), les arguments precedents ne sont pas applicables et 
on montre I’existence avec un argument d’approximation. 

Keywords: Mean Field Games, density constraints, variational formulation, convex 

duality methods 
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1. Introduction 

The theory of Mean Field Games (shortly MEG in the sequel) was introduced recently 
and simultaneously by J.-M. Lasry and P.-L. Lions ([2, 3, 4]) and M. Huang, R. P. Mal- 
hame and P. E. Gaines (see [5]). The main objective of the MEG theory is the study of 
the limit behavior of Nash equilibria for symmetric differential games with a very large 
number of “small” players. In its simplest form, as the number of players tends to inhn- 
ity, limits of Nash equilibria can be characterized in terms of the solution of the following 
coupled PDF system: 

-dtu{t,x) - uAu{t,x) + H{x,Vu{t,x)) = f[m{t)]{x) in(0,T]xM'^ 

^ dtm(t,x) - uAm(t,x) - div (VpH(x,Vu(t,x))m(t,x)) = 0 in(0,T]xM'^, 

m(0,a:)=mo, u{T,x)=g{x) in M'^, 

(MEG) 

where i7(x, •) is convex. The Hamilton-Jacobi-Bellman (HJB) equation in (MEG) charac¬ 
terizes the value function n[m] associated to a stochastic optimal control problem solved 
by a typical player whose cost function depends at each time t on the distribution m{t, ■) 
of the other agents. We remark that this interaction can be global, e.g. if /[m(t,-)](x) is 
a convolution of m(t, •) with another function, or local, i.e. when f[m{t)]{x) can be iden- 
tihed to a function f(x,m(t,x)). The Fokker-Planck equation (FP) in (MEG) describes 
the evolution m[u] of the initial distribution mo when all the agents follow the optimal 
feedback strategy computed by the typical agent. We refer the reader to the original 
papers [2, 3, 4] and the lectures [6] for more details on the modeling and the relation with 
the system (MEG). See also [7, 8] for a survey on the subject. 

For local couplings /(■,m), system (MEG) can be obtained (at least formally) as the 
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optimality condition of problem 


min 

s.t. 


dtm - uAm + div(w) = 0, m(0,x)=mo, 



dx, 


( 1 . 1 ) 

rm 

with F(x,m) := / /(x,m')dm', L(x,v) := H*(x,v) (where the Fenchel conjugate 

Jo 

H*{x,v) is calculated on the second variable of H) and niQ e satisfying that 

mo > 0 and / modx = 1. This type of approach, including also the degenerate first 

order case (z/ = 0), has been studied extensively in the last years in a series of papers 
[9, 10, 11, 12]. The optimization problem above recalls the so-called Benamou-Brenier 
formulation of the 2-Wasserstein distance between two probability measures, which gives 
a fluid mechanical or dynamical interpretation of the Monge-Kantorovich optimal trans¬ 
portation problem (see [1, 13]). We refer the reader to [14], [15] and the recent work 
[16] for some optimization methods to solve numerically (MFG) based on the formulation 
( 1 . 1 ). 

With a well-chosen time-averaging procedure, one can introduce stationary MFG sys¬ 
tems as an ergodic limit of time dependent ones (see [17, 18]), 


-zzAm(x) + H{x, V'u(x)) - A 


-z/Am(x) - div {'VpH(x, Vu(x))m(x)) 




0, m > 0. 


/(x,m(x)) in 
0 in 


(MFG^) 


At least formally (MFGoo) can be obtained as the first order optimality condition of the 
problem 


min 

s.t. 


w(x) 


/ ]m(x)L X,- 
Jud- ( y m(x) ^ 

- Am + div(tx) = 0, / m(x)dx 


+ F{x, m(x))| 

= 1 , 


dx, 
m > 0. 


( 1 . 2 ) 


The existence of smooth solutions for both evolutive and stationary MFG systems has 
been obtained in various settings in a series of papers (see for instance [19, 20, 21, 22] 
and the references therein). The used techniques combine variational arguments and 
sharp PDF estimates. Connections between stationary MFG systems and the so-called 
Evans-Aronsson problem have been also recently studied in [23]. 

The objective of this work is to rigorously study the optimization problem (1.2) 
with the additional constraint m < 1 a.e. Formally this should be linked to a system 
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like(MFGoo) with m <1 a.e. and an additional Lagrange multiplier corresponding to the 
new constraint. Moreover, in view of the interpretation of (MFG) as a continuous Nash 
equilibria, we expect that our derivation of an MFG system with a density constraint is 
linked to symmetric games with a large number of players on which “hard congestion” 
constraints are imposed. Similar models in the framework of crowd motion, tumor growth, 
etc. have been already studied in the literature (see for instance [24, 25]). In the case of 
MFG systems, we refer the reader to the papers [26] (for evolutive systems) and [21] (for 
stationary systems), in which “soft-congestion” effects, meaning that people slow down 
when they arrive to congested zones, are studied. Let us remark that in [6] it is also ex¬ 
plained how to study systems like (MFGoo) by means of a (degenerate) elliptic equation 
in space-time. However, this approach with the additional constraint m < 1 a.e. seems to 
be ineffective. 

The question of hard congestion effects/density constraints for MFG systems was hrst 
raised in [27]. More precisely, in the cited reference the author asks if a MFG system 
can be obtained with the additional constraint that the density of the population does 
not exceed a given threshold, for instance 1. To the best of our knowledge, this work is 
the hrst attempt to investigate this question. The stationary setting plays an important 
role in our study and we expect to extend our results to the dynamic case in some future 
research. 

Let 12 c (d > 2) be a non-empty bounded open set with smooth boundary and such 
that the Lebesgue measure of 12 is strictly greater than 1. Moreover, let / : 12 x M ^ M 
be a continuous function which is non-decreasing in the second variable and dehne iq : 
M X ^ M (with R := R u {+oo}) and Cg : hF^’'?(12) x ^ R as 




(1,3) 


+ 0 O, otherwise. 


We consider the problem 



where, as before, F{x,m) is an antiderivative of f{x,m) with respect to the second 
variable. We divide our main results in two classes, depending on the value of q. 

Case 1: q > d. In this case, using the classical direct method of the calculus of 
variations, we prove the existence of a solution {m,w) of {Pq)- Using that m e hF^’'^(12) ^ 
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C(fi), we are able to compute the subdifferential of Cr('m,w) for any 1 < r < g. It seems 
that this type of result is new in the literature. Moreover, the continuity of m allows us 
to prove that the constraints in (Pg) are qualified (see e.g. [28, Chapter 2]). Using the 
computation of the sub differential with r = q and classical arguments in convex analysis, 
we derive the existence of n € (s e [l,d/{d- 1)[), A e M and two nonnegative 

regular measures p and p such that 


-Au ++ p - p - X = f{x,m), 
-Am - div Vn] = 0, 


Vm • n = Vu • n = 0, 

J mda: = l, 0<m<l, 
n 

spt(p) £ {m = 0}, spt(p) £ (m = 1}, 


in 

in n, 
on dVL, 

in fl. 


{MFGg) 


where the system of PDEs is satisfied in the weak sense, ‘spt’ denotes the support of a 
measure and g' := g/(g - 1). In the above system, p appears as a Lagrange multiplier 
associated to the constraint m < 1 and can be interpreted as a sort of a “pressure” term. 
We also compute the dual problem associated to (Pg) recovering (MFGg) by duality. 
Finally, in the open set {0 < m < 1} we prove some local regularity results for the pair 
(m, u). 

Case 2: 1 < q < d. In this case, even if the existence of a solution still holds true, m is 
in general discontinuous, which implies that the arguments employed in the computation 
of the sub differential of Cq{m,w) are no longer valid. Moreover, the discontinuity of m 
implies that the constraint 0 < m < 1 is in general not qualified. In order to overcome 
these issues, we use an approximation argument. By adding the term eCr{m,w) with 
r > d to the cost function and using the arguments in Case 1 we obtain a system similar 
to {MFGq) depending on e. Then, by means of some uniform bounds with respect to e 
and recent results on estimates on the gradients for solutions of elliptic equations with 
measure data (see [29]), as e i 0 we can prove the existence of limit points satisfying 
(MFGq) where the concentration properties for p and p have to be understood in a weak 
sense. 

The structure of the paper is as follows: in Section 2 we set the basic notations 
and prove some preliminary results including the computation of the subdifferential of 
Cq{m,w). In Section 3 we define rigorously problem (Pg) for the case q> d and we prove 
the existence of a solution as well as the qualification property of the constraints. In 
Section 4 we characterize the solutions of (Pg) in terms of (MFGq) still in the case q> d. 
Moreover, we prove some local regularity results and we derive the dual problem. The 
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uniqueness of the solutions is also discussed. In Section 5 we complete the proof of the 
previous statements for any 1 < g < d by means of an approximation argument. Finally, 
in the Appendix we recall some important results about elliptic equations with irregular 
right hand sides. 

Acknowledgement: Both authors would like to thank F. Santambrogio for his con¬ 
stant interest and valuable discussions during the preparation of this paper. They are 
particularly thankful for his suggestion of the approximation procedure to treat the less 
regular cases. The hrst author would like to thank the hospitality of the University of 
Limoges several times during this project. He also thanks the hospitality of the Fields 
Institute, Toronto, in the framework of the Special Semester on Variational Problems in 
Economics, Physics and Geometry, Fall 2014. Both authors thank the partial support of 
the project iCODE: “Large-scale systems and Smart grids: distributed decision making” 
- “strategic crowds”. F. J. Silva benehted also from the support of the “FMJH Program 
Gaspard Monge in optimization and operation research”, and from the support to this 
program from EDF. Last but not least, we thank L. Brasco for the discussions on the 
Calderon-Zygmund theory, directing us towards the references on this topic used in this 
paper. 

2. Notations and Preliminary Resnlts 

We hrst hx some standard notation. Let H (d > 2) be a non-empty, bounded 
open set with a smooth boundary, satisfying a uniform interior ball condition, and denote 
by n the outward normal to dVt. Let us set | • | for the usual euclidean norm on and, 
given a Lebesgue measurable set A £ if it is not ambiguous, we also use |A| for its 
d-dimensional Lebesgue measure. 

We denote by the space of (signed) Radon measures dehned on H. We set 

(respectively 07t_(H)) for the subset of 9IT(H) of non-negative (respectively non¬ 
positive) Radon measures. Given the Hahn-Jordan decomposition m = -m~, with m'^, 

m~ € 9Jt+(H), we set := nPikl) + m~{VL) for the total variation of m. We also denote 

by 91lac(f^) and ^^^(H) the spaces of absolutely continuous and singular measures w.r.t. 
the Lebesgue measure, respectively. For notational convenience, if m e IHac(H) we will 
also denote by m its density w.r.t. the Lebesgue measure. Given p. e OJt(H) we set /iLA 
for its restriction to A £ H, dehned as n\—A(B) := fi{AnB) for all B e B{kl) (where R(f2) 
denotes the Borel cr-algebra on H). Finally, given A £ M.^, we set Xa(x) = 0 if a: e A and 
+00 otherwise. Moreover we set 1^(3^) = 1 if a: e A and 0 otherwise. 
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Now, let g > 1 be given and set g' := g/(g - 1). Consider the sets 

Aq! := {{a,b) e M X : a + < 0}, 

Aq'■= {{a,b) € X (a(x),b{x)) € Aq', for a.e. x € f2} , 

and recall the functions iq and Cq dehned in (1.3). We have the following result. 

Lemma 2.1. Suppose that q > d and let 1 < r < q. Then, the following assertions hold 
true: 

(i) The closure of Ar' in x L'^'(Tiy is given by 

Af^ = l^{a,^)€m(n)xu\ny a + ^i/^i^'^oj, ( 2 . 1 ) 


where the inequality in (2.1) means that for every non-negative (f € (7(12) we have that 

1 


f(j)(x)da(x) + — f (j)(x)\(3(x)\'^' dx <0. 
Jn r' Jn 


( 2 , 2 ) 


(ii) Restricted to W^’'J(f2) x the functional is convex and l.s.c. Moreover, for 

every {m,w) e iy^’^(12) x Li{Tiy, it holds that 


Cr{rn,w) = sup / [a{x)m{x)-i-ft{x) ■ w{x)] dx 

\ r r ] (^- 3 ) 

= sup _m{x)da{x)+ / (3{x)-w{x)dx 

(rvA^G A , JVl 


and C*{a,l3) = ^ (hh^’'J(12))’^ x L^'{VlY. 

Proof of (i). Let {an,l3n) ^ Ar' be a sequence converging to some (a,/?) in (VL^’'?(f2))* x 
L7(12)'^. Then, for any non-negative (f> e 1T^’'J(12) we have that 

0(x)a„(x)dx<— f (/)(x)|/3n(x)r dx. 

r' JQ 

Since in L^'(QY, except for some subsequence, |/5n(x)|^' ^ |/5(x)|^' for a.e. x € 12. 

Having positive integrands in the second integral, by Fatou’s lemma we obtain 

{a,(j))^w^,gy wi,g<-^J^(j){x)\l3{x)fdx for all 0 e ^^’^(ll), 0 > 0. 

In particular, letting 0 = 1, we have that jd € V' {VlY and by [30, Chapitre I, Theoreme V] 
we can extend a to a linear functional over (7(12), i.e. to an element in 911(12), satisfying 
(2.2). This proves one inclusion in (2.1). 
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In order to prove the converse inclusion, let (a,/?) be an element of the r.h.s. of (2.1). 
Equivalently, 

< 0 a.e. in Q and a® < 0 , 

where and a® denote the absolutely continuous and singular parts of a with respect 
to the Lebesgue measure, respectively. We shall construct different approximations for 
and (3 on the one hand and for a® on the other hand. For 7 > 0 and x e we set 
B^(x) = (y e ; |^ - x| < 7 }. Consider a molliher 7 : R*^ R satisfying that rj e C'^(R'^), 

7 > 0, r r]{x) dx = 1, spt( 7 ) c Si(0) and r]{x) = r]{-x) for all x e R'^. Now, for e > 0 set 

rig := {x e n : dist(x,912) > e}, 7 g(x) := —^(x/e), 
and for all x e 12 and i = 1 ,..., d, let us dehne 

dg(x) := / 7 g(x-?/)a“"(?/)dylQ^(x), ^^x) := / T]e{x - y)I3\y)dyln,{x). 

Jn Jn 

By convexity and Jensen’s inequality, for all x e f2g we have that 

^ * Ve{x) < 0, (2.4) 


and so (Q!g,/3g) € Ar' and one easily checks that dg 0 ;“^ in 9Jl(f2) and /3g ^ in L^'^QY. 

In order to approximate a® let us define the following kernel: for x € 12 and e > 0 let 
us set Pe {^B^{x)nn) x) n f2|. Note that for all x € 12 we have that dx in 0Jt(f2) 
as £ i 0. Given ?/ e 12 and £ > 0 let us define 

de(y)--= rpg(?/)dW(x). 

Jn 


Observe that for all £ > 0 the function dg is non-positive and, due to our regularity 
assumption on 912, we have that dg € L°°(f2). Let us show that dg ^ a® in 921(12). For 
any cj) e (9(12), Fubini’s theorem yields 


f <P{y)^e{y) dy = f (j){y) rp^(?/)da"(x)dy 

Jn Jn Jn 

_(j){x) da^(x) as £ i 0 , 

Jn 


n Jn 


0(d)P£(l/)d?/da"(x) 


where we have used that cj) is uniformly continuous in 12 (since this set is compact) and so 

/ (j){y)pAy) dy (j)(x) uniformly in 12 as £ i 0 . 

Jn 


This proves the convergence of dg. Dehning, we have that (ae,/?^) e Ar' and 

^ (a,/S) in 97t(r2) x The embedding 97t(r2) ^ (hh^’^(r2))* implies that 

the convergence also holds in (hh^’'?(r2))’" x from which assertion (i) follows. 

Proof of (ii). It suffices to show (2.3) (here we remark that by the Sobolev embedding we 
identify m with an element in (7(12), hence the second integral is meaningful). Indeed, 
(2.3) shows that Cr is the supremum of linear and continuous functionals, hence it is 
convex and l.s.c. For A: € M set Ar'^k '■= {{ci,b) ^ ; o k -k, maxj=i |6j| < k}. Since a.e. 

in 12 we have that 


lim sup {am(x) + b ■ w{x)} = sup {am(x) + b ■ w{x)} , 

(a,b)eA^i f. (a,b)eA^i 

and (a,b) = (0,0) € Ar^^k, by monotone convergence we have that 

£r.(m,tn) = lim / sup [am{x) + b ■ w{x)]dx. (2.5) 

Note that if \{m < 0}| > 0, then by (2.1), we readily check that both sides in (2.3) are 
equal to +oo. On the other hand, note that for every {m,w) € M x with m > 0 there 
exists a unique pair (a(m, rc), 5(m, tc)) e M x such that 


sup {am + b ■ w} = a{m,w)m + b(m,w) ■ w. 

(a, 6 )e^^, J, 


Indeed, 


- —m + b-w\ and a{m,w) = 


\b(m,w)\'^' 




( 2 . 6 ) 


which are well-dehned by the strict concavity of the objective function. Moreover this im¬ 
plies that M+xM'^ 3 (m, w) (a(m, ta), b(m, w)) € MxM'^ is continuous and measurable and 
thus hF^’^(12) X (L'?(12))'^ 3 (m,w) (a(m, tc), 5(m, re)) e L°°(12) xL°°(12)‘^ is well defined. 

Therefore, defining Ar',k ■= {{o,,b) € L°°(12) x L°°(12)'^, (a{x),b{x)) e Ar'^k for a.e. x e 12} 
we get that 


/ sup [am{x) + b ■ w{x)]dx 

(a,b)^A^, ,^ 


sup 

{a,b)^Ar 


[a(x)m(x) + 6(x) ■ tc(x)] dx. 


which, together with (2.5), implies that 

Cr{fn^w) = \im. sup / [a(x)m(x) + 5(x) • m;(x)] dx 

[a{x)m{x) + b(x) ■ w(x)]dx, 


= sup 

(a,b)eAr 
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proving the first equality in (2.3). The second equality follows from (i) and the continuity 
of the considered linear application. Finally, the identity C* = ^ is a consequence of (i) 

and (2.3). □ 


Remark 2.1. We refer the reader to [31, Chapter 5] for the proof of the semicontinuity 
of Cr in a more general setting. 

For m € denote := jx € hi: m{x) = 0| and = jx € hi: m(x) > 0|. Note 

that ii q> d, then E[f is closed. 

Theorem 2.2. Let (m,w) € x L‘i{VtY (q > d) and 1 < r < q. Suppose that 

Cr(m,w) < oo. Then, if v := {wjmflE^ i we have that dCr{m,w) = 0 . Otherwise, 

Cr is subdifferentiable at (m,w) and 

aC(m,zx) = |(a,/3) € and /3LEf^ = \v\^-%y (2.7) 


In particular, the singular part of a is concentrated in E^ ■ 

Proof. First note that since Cr{m,w) < oo, we have that \{m < 0}| = 0 and ta = 0 a.e. in 
E^f. By Lemma 2.1 for all {m,w) e x m > 0 we have that 


dC 


.(m,tc) =argmax(„^^),;^| f_mda+ f /3-wdx 


We claim that 


sup 

{a,g)LAC 



sup 

0eL^'(n)d- 



wdx. 


( 2 . 8 ) 


Indeed, the inequality “>” is immediate. To show the converse inequality for every e > 0 
let (as, fie) ^ such that 


_m(x)das{x)+ / fi,;{x) ■ w{x) dx > sup t _m{x)da{x)+ / fi{x) ■ w(x) dx\ - e. 
JQ Jfl {a,p)LAf' i 


Then, denoting by s the r.h.s. of (2.8), by (2.1) and the previous inequality we have that 

s >— [ m\fis\^'dx + f fi^-wdx> sup | fm{x)da{x)+ f ,d(x) • t(;(x) dx| - e 
r' Jn Jn (a,0)ap i 
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and so (2.8) follows by letting e ^ 0. Let us prove now that ii v i L^{E^Y, then 
dCr(m, w) = 0. We argue by contradiction supposing that there exists (a, /5) e dCr{m, w). 
By (2.8) and the assumption tc = 0 a.e. in E^, (3 must be a solution of the problem 


inf J(/3), where J(/5) := [ \-\PY'- v 13 

Jnlr' 


mdx. 


{ny 

Since vm = ta e q>r, we have that J is Frechet differentiable and 

O = DJ0)I3= [ r|/3r'-2/3-nl-/3mdx for all/5 e L"'(f2)^ 

JEy ^ 

which implies that, since (3 is arbitrary and m > 0 on 

v{x) = for a.e. x e F'™, 


(2,9) 


( 2 , 10 ) 


( 2 , 11 ) 


which is a contradiction because ‘^f3 e L'’(£'{")'^. Now, assume that v € L^^EY'Y 
let us prove that dCr{m,w) y 0. Define the functional J : L^(£'i)‘^ ^ M (where L'Y^iE^') 
denotes the space of measurable functions dehned in E^^ which are integrable w.r.t. the 
measure m) as 


J(/3) 



Since r' > 1, we have that J is coercive, continuous and strictly convex. Since LY^{EY^)^ 
is a reflexive Banach space, classical results in convex analysis imply the existence of 
a unique j3 e LY^^EY^Y such that J0) = inf{J(/5) : (3 e L^(i?{")'^}. The first order 
optimality condition implies that f3 satisfies (2.10)-(2.11) and so 


j3{x) = |n(x)|^ ‘^v{x) for a.e. x € E^- 


( 2 . 12 ) 


Since v € L'"{EYY have that (3 € L^^'^E^Y- Moreover, using that L^'{E^Y - 
relation (2.8) implies that 0'0) ^ dCr(m,w) and so dCr(m,w) y 0 . Now, let 

{a,j3) € dCr{m,w). The expression for Ar' in (2.1) implies that (-(l/r')|/5|^\ ^) attains 
the supremum on the r.h.s. of (2.8). Therefore, we must have 


[ mda + — [ \BY mdx = f mda + — f \BYmdx = 0. 
Jn r' Jn JEy r' JEy 


(2.13) 


Let us prove that a\—E'Y is absolutely continuous w.r.t. the Lebesgue measure restricted 
to EY- Let B e B{EY), such that \B\ = 0. Then, (2.13) implies that 


f mda + — f \f3Y'mdx+ [ mda 
JEV^^B r' JEV^^B JB 


0 . 
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By a standard argument using Lusin’s theorem (to approximate the by continuous 

functions) and ( 2 . 1 ) we must have that f^mda = 0 and since m > 0 on i?™ we conclude 
that a{B) = 0. Thus a\—E^ « | •|L£'™. In particular, spt(a^) £ and, denoting still by 
a the density of a restricted to i?™, a{x) + (l/r')|^(x)|^' < 0 for a.e. x e E'^. Therefore, 
by (2.13) we have that 



dx = 0, 


and since m > 0 on we conclude that a = a.e. in E'^. Using (2.8) we get that 

jS solves problem (2.9) and so f3 = a.e. in E^ from which the result follows. □ 


Remark 2.2. Note that redefining the domain of Cr as C{Q) x L'?(fl)‘^, the above proof 
shows that the conclusions of the Theorem 2.2 are still valid in this setting. 


3. The optimization problem 

In this entire section we suppose that q> d. Let / :flxR^Mbea continuous function 
in both variables and increasing in the second variable. Let us dehne the function 

nm 

X M 3 (x, m) 1 -^ F(a:, m) := / f{x,s)dseM. 

Jo 

Note that for every hxed x e the function m F{x,m) is convex. Let us dehne 


.F : lU^’^(fl) ^ M as F’(m) := [ E{x,m{x))dx. 

Jn 

Given w € we consider the following elliptic PDE 

( -Am + div(t(;) = 0 in G, 

( {ym-w)-n = 0 on clf 2 . 

We say that m e W^’^iVL) is a weak solution of (3.2) if 

/ Vm(x) ■ V</?(x)dx = / w{x)-V(p{x)dx V(/)€G^(G) 

Jn Jn 


(3.1) 


(3,2) 


(3,3) 


By Lemma 5.3 and Lemma 5.4 in the Appendix for a given w e L'?(f2)'^ equation (3.2) has 

a unique solution m € satisfying that / m(x)dx = 1. We consider the following 

Jn 

optimization problem; 


inf Jg(m,w) := Cq(m,w) + E{m), 
{m,w)iKp 


(A 
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where the set of constraints )Cp is dehned as 

/Cp := |(m,'u;) e x ; such that (m,w) satisfies (3.2), L m{x) dx = 1 , m < 1 


Remark 3.1. Since Cq = Cq + X{m>o}, the constraint m > 0 is implicitly imposed in {Pq). 

Given s e [ 1 , oo[ we set hh*’^(G) := jn e ■ J' = o|- Now, let us define A : 

^ and B : L<i(Qy ^ as 

{{Am, 0)) := f vm{x) • V0(x) dx, {{Bw, 0)) := - [ w{x) • V0(x) dx, 

Jn Jn 

for all m e w e and 0 e (G), where we used ((•,•)) to denote the 

duality product between (w^'^ (^)) and (G). Since A and B are linear bounded 

operators, the adjoint operators A* : (G) and B* : (G) -> L‘^'{iiy 

are well-dehned and given by 

{A*(j),m) = / V0(x) ■ Vm(x) dx, {B*(j),w)qf^q = - / V0(x) ■ tc(x) dx, 

w/ 

where we have used (•,■) to denote the duality product between and 

and {■r)q',q to denote the duality product between L^'{fiy and Li{flY. Now, let / : 

^ C{Q) be the Sobolev injection, which is well-dehned since q > d (see [32]), and 

let C:={ze C(n) ] z<l}. Let us set X := x L^ny, Y := (fyi’^'(G))" x M x G(G) 

and dehne the application G '■ X as 

G{m,w) := ^Am + Bw,J^ m(x) dx - 1, Jmj . 

By setting tC := {0} x {0} xC ^Y we have that /Cp can be rewritten as 

JCp = {{m,w)^X :G{m,w)^X}. 


Since A, B and / are linear bounded operators, we have that /Cp is a closed and convex 
subset of W^’^{VL) x L'i{VLy. 

Theorem 3.1. Problem {Pq) has (at least) one solution {m,w). 
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Proof. Since |r2| > 1 we have that {m,w) := (l/|r2|,0) belongs to ICp and the cost function 
is finite. Now, let {mk,Wk) e K,p be a minimizing sequence. Since m < 1 a.e. in hi 
and Cq{mk,Wk) is bounded uniformly in k, we get that ItCfcllLi is bounded. Therefore, 
there exists w e such that, except for some subsequence, Wk ^ w weakly in 

L^{QY. In addition, by Lemma 5.4 and the boundedness of Wk in L^{VlY we have 

that Vrrik is uniformly bounded in L^^Y- Since f mk{x)dx = 1, Poincare’s inequality 

Jn 

< C\\Vmk\\L<i) implies that is bounded in Thus, there exists 

m e W^''^{YL) such that, except for some subsequence, rrik m weakly in Using 

these convergences, we get that Am + Bw = 0 and / m(x) dx = 1. The continuous 

Jn _ 

embedding / preserves the weak convergence and C is weakly closed in C'(r2). Thus, 
Im e C, which implies that {m,w) € fCp. Since Jq is convex and l.s.c. w.r.t the weak 
topology in IU^’'?(r2) x L^(flY (by Lemma 2.1) we get that Jq{m,w) = inf{j7g(mi,tci) : 
(mi,rci) e )Cp}. □ 



Now, we prove a constraint qualification result for problem (Pg) (see e.g. [28, Chapter 
2]), which is crucial for deriving optimality conditions. We set dom(i7g) := {{m.,w) € 
X L^{VtY '■ Jq{m,w) < oo}. 


Lemma 3.2. We have that 


0 e int {G(dom(j7g)) - JC} . 


(3.4) 


Proof. We need to prove that for any given (hi,h 2 ,h 3 ) e Y small enough there exists 
(m,w,c) € dom(57g) xP such that 


Am + Bw = 5i, / m(x)dx = 1 + (52, I(m)-c = S 3 . 

Jn 


(S) 


We observe that (m,0) € dom(j7q), for all m € n dom(P) non-negative, which 

implies that we can search the solution of (S) in the form (m,0,c) e dom.{Jq) x C. First 
of all, note that for mo := l/|f2| we have that 

v4mo = 0, / mo(x)dx = l, Jmo = mo e int(C). 

Jn 

By Lemma 5.4 (see the Appendix), there exists mi € such that 

Ami = 5i and f mi(x) dx = 1 + ^ 2 . (3.5) 

Jn 
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Setting Sm := mo - mi, we obtain that A 6 m = -5i and / Sm(x)dx = - 62 - By Lemma 

Jn 

5.4 the linear bonnded operator b m ^Am, J m(a;)dxj e (^)) x M is 

surjective and so, by the Open Mapping Theorem, there exists h e VL^’'?(r2) such that 

Ah = 0 , J^h{x)dx = 0 and = O ^|| (^i, • 

In particular, a.s q> d the Sobolev inequality implies that 

||/(fe)-/Om)IU- =o(||W,iOII(„,y(a,)•,.). (3.6) 

Now, let r := h-5m and for 7 > 0 let us dehne ■= mi+jh, which by construction solves 
(3.5). Since mo e int(C), if 7 is near to one (and 61,62 are small enough) then, by (3.6), 
I{m^) = /(mi) + 'jl{ 6 m) + jl^r) € int(C). Thus, if ||(53 ||l~ is small enough we have that 
c := I{m^) - 63 e int(C). Thus, (m.y,0,c) solves (S) and e dom(j7g). The result 

follows. □ 


4. Optimality conditions and characterization of the solutions 

The purpose of this section is to derive optimality conditions for problem [Pq) and, 
as a consequence, to obtain the existence of solutions of (MFGq). As in the previous 
section we will assume in all the statements that q > d. Our strategy relies on a “direct 
method”, which uses the constraint qualihcation condition established in Lemma 3.2 
and the characterization of the sub differential of Cq (see Theorem 2.2). The uniqueness 
and local regularity of the solutions are also discussed. Moreover, in Subsection 4.1 
we formulate the associated dual problem, and we provide an alternative (but related) 
argument to derive optimality conditions. 

Let us dehne the Lagrangian H x L^{yLY x (O) x 911(0) x M ^ R as 


61{m,w,u,p, X) := J'q{m,w) - (lAm + Bw,u} + J Jmdp + A ^y'mdx-lj (4.1) 
= Jqim, w) - {A*u - I*p - X,m) - J B*u ■ wdx - A. 


Remark 4.1. Since the inclusion 1T^’'^(0) ^ 67(0) is dense, for every measure p € 971(0) 
the adjoint of the injection operator I*p at p can be identified uniquely with the restriction 
of p to 1T^’'?(0). Thus, for notational convenience we will still write p for I*p. 
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Recall that for a Banach space X and a convex closed snbset K ^ X, the normal cone 
to K at a: e K is dehned as 

Nk{x) ■■= {x’" € X* : {x*,z - x)x*,x < 0, Vz e K}, (4.2) 

where {■,-)x*,x denotes the duality pairing of X* and X. Using [28, Example 2.63] we 
have 

Nic{G{m,w)) = |(m, A,p) e (fi) x M x ; spt(p) n = 1 }|. ( 4 . 3 ) 

Now, we provide the first order optimality conditions associated to a solution (m, w) of 

(A)- 

Theorem 4.1. Let {m,w) € JCp be a solution of problem {Pg). Then, v := {w/m)l{rn>o} ^ 
L^iyiY and there exists {u,p,X) e lUl®(r 2 ) x x M for all s e]l,d/{d- 1)[ and 

(a,/?) 6 dCq{m,w) such that A*u e 911(12) and the following optimality conditions hold 
true 

' a + f(-,m) - A*u + p + X = 0, 

[3 = B*u, 

Am + Bw = 0, 

spt(p) £ {m = 1 }, 0 < m < 1 , / mdx = l, / ndx = 0 , 

JQ. Jn 

where the first equality holds m 911(12). Conversely, if there exists {a, {3,u,p, X) e dCg{m,w)x 
IU^’'J'( 12 ) X 91t+(12) xM such that (4.4) holds true, then (m,w) solves {Pq)- 

Proof. By Lemma 3.2 the problem is qualihed (see e.g. [28, Chapter 2]). Thus, by classical 
results in convex analysis (see e.g. [28, Theorem 2.158 and Theorem 2.165]), we have the 
existence of (n,p. A) e Nx.{G{m,w)) such that 

(0,0) € a(^^^)£(m,n;,M,p,A). (4.5) 

Since Cq is hnite at (1/|12|, 0) and the other terms appearing in £ are differentiable, by [33, 
Chapter 1, Proposition 5.6] and (4.5), we must have that Cq is subdifferentiable at {m,w). 
Thus, by Theorem 2.2 and (4.2) we get that v e there exists (a,/5) e Aq' such 

that (4.4) holds true, with the first equation being an equality in (1U^’'?(12))’". Since, 
except by A*u, all the other terms can be identified with elements in 911(12), we have that 
A*u can be identified to an element of 911(12). Using classical elliptic regularity theory (see 
[34, Theoreme 9.1]) we get that u € Wo’^(12) for any s e]l,d/{d- 1)[. The fact that (4.4) 
is a sufficient condition follows also by the convexity of the problem (see [28, Theorem 
2.158]). □ 


16 




As a corollary we immediately obtain the following existence result for a MFG type 
system with density constraints 

Corollary 4.2. There exists {m,u, iJ,,p, X) e x Wo^’^(r 2 ) x 9Jt+(r2) x 9 Jt+(r 2 ) x M 

(s e]l,d/((i - 1)[) such that 


' -Au + ^\Vu\i'+ p - p - X 
Vu ■ n 


-Am - div ^m\Vu\i-^ 

Vm • n 

/ mda: = 1, 

JQ 

spt(/i) £ {m = 0 }, 


f{x,m) in n, 

0 on d^l, 

0 in n, 

0 on d^l, 

0 < m < 1, in n, 
spt(p) £ {m = 1 }, 


{MFG,) 


where the coupled system for {u, m) is satisfied in the following weak sense: for all p e 

c^(n) 

f Vu-S/pdx + f —IVul^'pdx - X f pdx+ fpd(p-p)= f f(x,m(x))p(x)dx, 

Jn Jn r Jn Jn Jn 

Vm + m|Vti|~ • Vtpdx = 0. 


Li' 


(4.6) 


Let us define := {x € hi: 0 < m(x) < l}. Note that by the continuity of m, is 
an open set. 

Remark 4.2 (The uniqueness of the solutions). Assuming that the coupling f is strictly 
increasing in its second variable, the objective functional in (Pg) becomes strictly convex 
in the m variable (and the set Kp is convex). Thus, the function m e iy^’5(r2) in {MEGg) 
is unigue, which implies also the unigueness of w ^ In particular, Vu € 

is also unigue on Ef^- The first identity in (4.6) with p e Cl{E'!fi) implies the unigueness 
o/ A e M. If p ^ G){E'(^) we obtain 

f 'Vu-Vpdx + f —\'^uL'pdx - X f pdx- f pdp= f f{x,m{x))p{x)dx, 

Jn Jn r' Jn Jn Jn 

which together with the condition spt(p) £ {m = 1} yields the unigueness of p. Using [35, 
Theorem 3-4] we obtain that on E'fi u € {fil) is unigue up to an additive constant 
which may differ on each connected component of the set Ef^. In general, we cannot 
expect unigueness for p. 
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Now, let us discuss some interior regularity properties for the solutions on the open 
set i?™. Our approach is based on a bootstrapping argument. 

Proposition 4.3. There exists 70 £]0,1[ such that 

and (4.7) 

If in addition, f € x M) for j € {0,1} and some 71 e]0,1[, we have that for some 

72e]0,l[ 

(4.8) 

Proof. Step 1. We show that there exists k> d such that u € By the classical 

Sobolev embeddings, this implies that u e Cyff{E^) (for some 7 e]0,1[). Let ri e]q',d/{d- 
1)[. Since u e hL«’^^(r2) we have that \Vu\T e (Vt). The continuity of / and the density 
constraint on m imply that /(■,m(-)) e L°°{Q). Thus, denoting by := ri/q', classical 
regularity theory (see [36]) yields u e {Eff). In particular, the Sobolev inequality 

^ dS^ 

(see e.g. [32]) yields u e Wio’/ and so |V'n|'^' e . We easily check that 

82 '■= d5\lq'{d- hi) > hi and so n e W’ioc^(-E'™)- Let us define the sequence hj+i := 

Since hj+i - 6 i > {q' + d - dq')l(d - 6 i)q' and q' + d - dq' > 0 , after a hnite number of steps 
we get the existence of i’" > 2 such that k := h** > d and u e Wyf^{E^). 

Step 2. Let us prove that m e Cyff°{Eff) for some 70 e]0,1[. Since m € W^’'?(r2) and q> d, 
we already have that m is Holder continuous. Having u e this implies that 

Vn e C'loc hence m|Vn|^VM e for some 7 e]0,l[. Using a Schauder- 

type estimate (see [37, Theorem 5.19]) we get that m € {El^) for some 7 ' e]0,1[. 

Step 3. Using the above regularity for m, if / e Uj°’[]'^(r 2 xM), the local Holder regularity 
for and [36, Corollary 6.9] imply that u e {Eff) for some 7 " eJO, 1[. Finally, if 

/ e ClYY (f^ X K), the local Hblder regularity of Vm and of D‘^u imply that u e {Eff) 
for some 7 "' 

□ 


4-.1. The dual problem 

In order to write explicitly the dual problem we will need the following Lemma con¬ 
cerning the Legendre-Fenchel transform of E. 

Lemma 4.4. Let E be defined by (3.1). Then its Legendre-Fenchel transform E* '■ 
(VF^’^(r2))* M is given by 

[ +CX), otherwise. 
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where F* denotes the Legendre-Fenchel transform of F w.r.t. the second variable. 

Proof. The result is a consequence of [38, Section 2], □ 

We recall that given a Banach space X and a convex closed set K £ X, the support 
function ax ■ X* ^ M is dehned as 

•= sup(a;*', 'ix*^X*. 

x^K 


Proposition 4.5. The dual problem of (Pq) (in the sense of convex analysis) has at least 
one solution and can be written as 

min I f F*{x, a) dx + X+p(Q,)\ (PDq) 

(u,p,\,a)^Ku l-/o J 


where 


Fd ■■= ] (u,p, A, a) € Wy (fi) X X M X : A*u + - p - A < a 

I q' 

where the inequality has to be understood in the sense of measures. 

Proof. The dual problem of (Pg) can be written as 


max 

(tt,p,A)€ 

vrh'^'(o)xa)t(o)x][ 


inf £,{m,w,u,p, X) - aic(u, X,p) 


(4.10) 


where il is dehned in (4.1) and we recall that K. ■= {0} x {0} x C. The fact that we have a 
max instead of a sup in (4.10) is justihed by Lemma 3.2 and [28, Theorem 2.165]. Now, 

note that _ _ 

p{Ll) if p € 91l_^(r2), 


(Tk{u,X,p) = (Tc{p) = 
On the other hand, we have that 


+ 0 O otherwise. 


(4.11) 


inf £(m, tc, M,p, A) = - sup -2{m,w,u,p, X), 

("I'Hi) (m,w) 


= -sup l{A*u-p-X,m) + / B*u ■ wdx - Jq{m,w)\ - X, 

(m.w) V -I 


(m,w) 

= -j;(A*u-p-X,B*u)-X 
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Since there exists {m,w) € dom(£q) at which T is continuous (take for example {m,w) = 
(l/|r2|,0)), for any {a,j3) e x L'^'(iiy we have that (see e.g. [39, Theorem 

9.4.1]) 

= - a,/3) + J*^"(a)} , 

= infa.(wi.^(o))‘ {x^(« - /5) + 

= : a + -\py <aY 

where we have used Lemma 2.1 and Lemma 4.4. Let us prove that the above minimization 
problem has a solution. First, by (4.9) the integral functional is l.s.c. with respect to 
the weak-* topology of measures. Let us take a minimizing sequence an e L^{Q). There 
exists a constant C > 0 such that 

C> f F*{x,an{x))dx> f [an{x)y{x) - F{x,y{x))]dx, VyeL°°(fi). 

Jfi Jn 


By choosing y{x) = sgn(a„(x)) (which is equal to 1 if an(x) > 0 and -1 if not), we obtain 
that an is bounded in L^(r 2 ). 

Therefore, when the sequence is identified to a sequence of measures, we get a 
weakly-* convergent subsequence to some a e 971(12). The constraint is convex and closed 
with respect to this convergence, so by the lower semicontinuity of the objective functional 
we have that a is a solution and, by Lemma 4.4, a e 97lac(12) as well. Using this result 
and (4.10), (4.11) and (4.12), the conclusion follows. □ 

Using the dual problem, let us provide an alternative, but related way, to obtain first 
order necessary and sufficient optimality conditions. By Theorem 3.1 and Proposition 4.5 
we know that that there exist (m,w) € JCp and (u,p,X,a) € /Cp optimizers for (Pg) and 
{PDq) respectively. Moreover, since Lemma 3.2 implies that problem (Pg) is qualified, by 
[28, Theorem 2.165] problem {PDg) has the same value as problem {Pg). Therefore, 


Cq{m,w) + F{m) = - f F*{■,a) dx - X - ac{p), 

r 1 , (4.13) 

Am + Bw = 0, m < 1, / mdx = l, A*u - p - X - a +—\B*u\‘^ < 0. 

Jn a' 
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Using the above relations, we obtain 

0 = jF(m) + y' F*{-,a)dx + Cq{m,'w) + - p - X - a, B*u) + X + ac{p), 

= F(m) + f F*(-,a) dx + Cq(m,w) + C*(A*u - p - X - a, B*u) + X +ac(p), 

Jn 

+ ^ f mdx + acip), 

>- - P - ^ ^ ^ a ^ mCx . Mp), 

= {A*U,m) + {B*U,w)qrq, 

= {{Am, u)) + {{Bw, u)) = 0. 

This means that all the inequalities in the previous list are actually equalities. Thus, 

(i) F{m) + F*{a) = {a,m) and so, using the fact that F is differentiable on lU^’'?(f2), 
we have a = /(■, m). 

(ii) Cq{m, w) + C*q {A*u - p- X - a, B*u) = {A*u - p - X - a,m) + {B*u, w)qi^q, namely 

{A*u - p - X - a,B*u) € dCq{m, w) 

(hi) ac{p) = (p,™-)ot(o),c(o)> which implies that p e Nc{m). 

Using (4.3), (4.13) and (i)-(iii) we recover system (4.4). 

5. Treating less regular cases via an approximation argument 

In this section we provide the proof of the existence of a solution of a suitable form of 
(MFGq) when 1 < q <d. Note that given w € L'?(f2)'^ the solution m of (3.2) is in general 
discontinuous. Because of the constraint 0 < m < 1, this implies that problem {Pq) is in 
general not qualihed (see [28, Chapter 2]) and thus the arguments in the previous section 
are no longer valid. To handle this issue, we propose an approach which is based on a 
regularization argument. 

Let us £x 1 < g < d and r > d. For e > 0 dehne Jg^e ■ VF^’^(fl) x L'’(f2)'^ ^ M as 

Jq^e{m,w) := Jq{m,w) +eCr{m,w). 
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Following the arguments in the proof of Theorem 3.1, problem 


inf Jq^e{'m,w) 




(n,s) 


admits at least one solution (me,We). Since e problem {Pg^e) is qualihed. 

Moreover, since both Cq and Cr are continuous at (rh,w) := (l/|r2|,0), by [33, Chapter 1, 
Proposition 5.6] we have that 


d{Cq(m, w) + eCrim, w)) 


dCq{m, w) + edCr{rn, w) for all (m, w) € x P{VlY. 


Therefore, exactly as in the proof of Theorem 4.1, if we define Vg, := {wg.lmY'^E'^^ we have 
that Ve e L’^iytY and there exist (uegPcXe) e VFi’*(r2) x 0Jl+(r2) x M (s e]l,d/{d- 1)[), 
(c(e,q,/3e,g) £ dCq{me,We) and (a£,r,/3e,r) ^ dCr{me,We) such that 


a^^q + eae,r - A*U^ + f(x, mY) + Pe + 

Ye,q S" ^Ye,r 

Arrip + Bwp 


0 , 

B*u,, 

0 , (5-1) 

/ m£da; = l, 0 < < 1, sptfpg) £ (mg = 1}. 

Jn 


Now, for e > 0, let us dehne Fq^^, Gq,ei Hq^e • 


as 


FqAz)-.= -\z\^ + ^z\Y GqYz)--=^\z\^ + ^\z\Y and Hq^z) := Gq^VPlYz)). 


For notational convenience, we set Hq := Hq Q. Elementary arguments in convex analysis 
show that Hq^^ Hq uniformly over compact sets. System (5.1) can be written in the 
following alternative form: 

Proposition 5.1. There exists e 91l_(r2) such that 


S" Hq^Y Vrig,) oc^ Ag 

-Arrie + div (meVFq*e(-VWe)) 
Vm^ ■ n = 0 

0 < me < 1 , 

spt(pe) c {m^ = 1}, 


f{x,mY, in fl, 

0, in n, 

VMe • n = 0, on dfl. 

/ m^dx = 1 , 

Jn 

spt(d:£) g {rrie = 0}. 


(MFGqY 
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in 


Proof. By Theorem 2.2 we have that 

= = = KrV, = \v,r\. 

On the other hand, since Vm^ e we have that := S/Ff^{-S/u^) e L'^{VlY. Using 

that SiFqYve) = (^e,q + ^(^e,r = “V^e in Ff^ and that VT’g’g = ^Ffe, we get that Ve = 

F^P Therefore, there exists fe ^ ^^^h that spt(^£) £ EY" and a.e. in O 

l3e,q = \Vs\'^~‘^Ve + L and /3^^r = ^{^FqYVe) - ^e^q) = \Ve\''~‘^Ve - Y/£)fe- 
Using the convexity of ^| ■ Y and ^| ■ |^', we easily check that 

and ^\I3e,r\’' ^ eY - Ve ’ L■ 


Hence 




q r q 

with an eqnality a.e. in In particular, we have the existence of a positive measure 
7 £ such that spt( 7 £) £ spt(^e) £ FY" and 

-Y>.P - ^\M' = - le. 

q r 

Since the dehnition of (ae,q,PE,q) and {ae,r,P6,r) implies the existence of two positive 
measures ae,q and aE,r such that spt(d£,g) £ FYP spt(Q;£,r) £ FY‘' and 


a, 


1 , 1 , 
£,q ~ ~ ^s,q and Cle^r — ~ 


Of 




the result follows by setting '■= -&e,q - ^^£,r - 7e- 
Now we present the main theorem of this section. 


□ 


Theorem 5.2. There exists {m,u,p, fi, X) e 1U^’'?(0) x kUi’^ (O) x 0Jt+(r2) x 0Jt+(r2) 
such that 

-Am + + p - p - A = f{x,m) in O, 

Vm • n = 0 on dfl, 

, 2-9 


-Am - div = 0 in O, 

Vm • n = 0 on dQ, 

J mdx = l, 0 < m < 1, in O, 

o 


{MFGq 
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where the coupled system for {u,m) is satisfied in the weak sense {see (4.6)). Moreover, 
defining 






mdx - 


f 

Jn 


\/m ■ S/udx 


(5.2) 


we have the inequality 



p, m) < 0. 


(5.3) 


Proof. Step 1: Bounds for A^, p^ and Note first that the second equation in (5.1) 
and Theorem 2.2 imply that a.e. in Also, in the set we have that 

^Fq,e{vfi) = -VWe and so in the identities v,. = VFffi-Vufi) and Hq^^{-\/ufi) = Gq^^{vfi) 
hold true. Now, by the second and third equations in (5.1) we get that 


/ Vu^-Vm^dx= / VUi.-Wedx = - / Fq ^{vfij-Vi^dx 

JQ Jn Jn 



‘^+e\vfif)dx (5.4) 


and so taking rUg as test function in the hrst equation of (5.1), we obtain 

\e+ I msdpe= / {Gq^e{ve)'me + VUe-Vm^- f{x,mfi)mfi) dx 

Jn Jn 

= / 1 — e-\veYme-f{x,mfi)mA dx, 
Jn\ q r / 

which implies that 


Xe+ medpe =-Cq{me,w^) - eCr{ms,w^) - / f{x,me)medx. (5.5) 

Jn Jn 

The optimality of {'ms,w^) yields 

0 < Cq{me,We) + eCr{m,,We) < Jq^e{^/\n\,0)-F{m^). (5.6) 


Thus, since / is continuous, 0 < rrig < 1, (5.5) and the fact that spt(p£) £ {m^ = 1} yield 
the existence of a constant Ci > 0 (independent of e) such that 

- Cl < Ae + / dps = Xe + IPsItv < Cl- (5.7) 

Jn 

On the other hand, by taking l-m^ as test function in the hrst equation of (5.1), a similar 
computation using (5.4) yields 

(|0|-l)Ae-|d£|Ty = / Hg^e{-VUe)dx + Cq{me,We)+sCr{me,We)- /(x,me,)(1 -m^)dx, 

w/ 

(5.8) 
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from which 


(|f]| - 1)A£ - Ittelry > - / /(x,me;)(l-mjdx > C2, 

Jn 


(5.9) 


where C 2 > 0 is independent of e. Since |r2| > 1, inequalities (5.7)-(5.9) imply that Ag is 
uniformly bounded w.r.t. e and so and are uniformly bounded w.r.t. e in 971(17). 

Step 2: Convergence ofVUe and 171^. By (5.8), as a function of £ we have that //^^^(-Vne) 
is uniformly bounded in L^{Q) which implies that is bounded in and that 

-Atig is bounded in 971(17). On the one hand, the boundedness of Ue in hh^’'J'(17) implies 
the existence of u e iy^’'?'(17) such that up to some subsequence converges weakly to u 

in Vhi’ 9 '(( 7 ). In particular, f udx = 0. On the other hand, the boundedness of -Au, in 

971(17) and [29, Theorem 1.3 with p = 2] imply the existence of s e]0,l[ and Aq > 0 such 
that VMe is uniformly bounded in lTj*’j[’^^‘’(17)‘^. By [40, Corollary 7.2] we can extract a 
subsequence such that VWe ^ Vm a.e. in 17 and so Hq^^{-S/Us) ->■ a.e. in 17. 

Now, in order to establish the convergence for iris, note that inequality (5.6) and the 
fact that 0 < rUg < 1 imply that Ws is uniformly bounded in L^(17)'^ for all e > 0. This 
means that, up to some subsequence, Ws is converging weakly in L'?(17)'^. Since Lemma 
5.4 implies that || Vmg||x ,9 < C'||tCg||L<j (for a constant C > 0 independent of e), by Poincare’s 
inequality we get that rUg is uniformly bounded in 1T^’'^(17). Extracting a subsequence 
again, there exists m such that rUg converges weakly to m in VP^’'?(17). By the compact 
Sobolev embedding, we get strong convergence in which implies that 0 < m < 1 

a.e. in 17 and / m dx = 1. 

Jn 

Step 3: The limit equations. The weak formulation of the second equation in [MFGq^s) 
yields 

/ Vmg-V(/?dx = - / msWF* Us) ■ Vpdx, for all 9 ? e C'°°(17). 

Jn Jn 

Since, extracting a subsequence, Ws = msVFj{-VUs) converges weakly in to some 

w, the weak convergence of rUg in 1T^’'?(17) implies that 

/ Vm-V(/?dx= / w-Vipdx, for all (^ € (7°°(17). 

Jn Jn 

Moreover, extracting a subsequence again, we get that 

2-q 

mg(x)vT]j(g(-VMe(x)) ^ -m(x)|Vn(x)|'j-i Vm(x) for almost every x e 17. 


,2-9 


The latter equality and Egorov’s theorem imply that w = Vtt from which the 

second equation in [MFGq) follows. 
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On the other hand, the weak formulation of the first equation in {MFGg^^) reads 

/ VUe-S/(pdx+ / Hq^^(-VUe)(pdx-X^ / ifdx- / ifd{pe + de)= / f(x,ms{x))(p{x)dx, 
JU JQ JQ Jo, JQ 

( 5 . 10 ) 

for any test function p € (7^(0). The continuity of / and the dominated convergence 
theorem imply that 

lim f f{x,me(x))ip{x)dx= f f{x,m{x)Mx)dx for all (p € 0 ( 0 ). 

JQ Jn 

The previous steps imply that we only need to study the limit behavior of the second 
term in (5.10). Since Hg is bounded in L^(VL), there exists 7 € 911(0) such that, 

extracting a subsequence, for all p € 0(0), f Hg ,{-Vu,)ipdx ^ T (pdy. Patou’s lemma 
implies that 

1 


f —\Vu\'^'(fdx f Hg ^(-\/Us)(pdx = f (pdy V (peO(O), yj > 0. 

Jn q' e^o Jn ’ Jn 

Defining, p € 911+(O) as dp := dy - ^Iv^l'^'da:, we obtain that 

f (pdp+ f —\Vu\‘^ pdx = \im f HqJ-Vu£)pdx for all (p € (7(0). 

Jn Jn a' e^o Jn 


( 5 , 11 ) 


Thus passing to the limit in (5.10) as e 0 we get 
1 

q' 


[ Vu-V(pdx+ r —|Vm|^ <pdx - A f pdx- f pd{p + d-p)= f f{x,m{x))p{x)dx. 
Jn Jn q' Jn Jn Jn 


Setting, p := p - d e 911+(O) we obtain the weak form of the first equation in [MFGg). 
Step 4- Proof of (5.3). By (MFGg^^) and (5.4) we have 

0 = / {l-m£)dp£- / m£dd£= / dp^ + / m^di-de-pe), 

Jn Jn Jn Jn 

= dp£+ / [Xe + f{x,me) - Hg^£{-Vu£)]m£dx - / VUe-Vrriedx, 

Jn Jn Jn 

rUg dx. 


f dpe+ f 

Jn Jn 

f dpe+ f 

Jn Jn 


+ -|xg|'’ + Ag + /(x, rUg) 

iq r 

1| 

Lg 


'^£1'' +Ag + /(x,mg) 


TTT-g dx. 
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By Fatou’s lemma we have 


/"-IVwl'^ mdx < liminf [ 
Jn a e^o Jn 


dx. 


Thus, letting £ ^ 0 and using (5.2), we get 

0> f dp + /" [-|Vn|'^'+ A +/(x,m) mdx 
Jn Jniq 

= / dp+ / VM-Vmdx+ / |Vm|'^ mdx + (/i-p, m). 

JQ, Jo, 

By taking u e (12) as test function in the second equation of {MFGq) we obtain that 
J Vm • Vn + m|Vn|~|Vnp dx = [vm • Vn + m|Vti|'^'] dx = 0, 


from which (5.3) follows. □ 

Remark 5.1. Inequality (5.3) is a sort of “weak concentration property”. In fact, by an 
approximation argument it is easy to see that we can take C{Il) n hF^’'?(f2) for the space 
of test functions in the first equation of (MFGq). Thus, if m is continuous, we would 

have that {p - p,m) = / md{p-p) and so (5.3) would imply that 

Jn 



0 and 



m)dp = 0, i.e. spt(/i) £ {m = 0 } and spt(p) £ {m = 1 }, 


as in Corollary 4-2. 


Appendix 

In this section we recall some classical results about the regularity of solutions of 
elliptic equations with irregular r.h.s. Recall that we set ((•,■)) for the duality product 
between {q > 1) and (12). The following surjectivity result holds true. 

Lemma 5.3. For any f € {Wl'^ (^))’" weak formulation of 

div(F) = / in 12, F • 77 , = 0 in 512, i.e. - f F(x) ■ Vpix) dx = (if,p)) (5.12) 

Jn 

for all ip € (12), has at least one solution F e Li{Tiy. 


27 





Proof. Let us consider the problem 


f IVwl'^'da:-((/,«)). 

1 . 9 ' o' Jn 


mm 
ueWp^' *? 


Since the cost function is strictly convex, coercive and weakly lower semicontinuous, we 
have the existence of a unique u e (hi) such that 

f \Vu(x)\^'~'^Vu{x) ■ V(j){x) dx = ((/, (p)) V (p e (hi). 

Jn 

The result follows by defining F = -\Vu\‘^'~^Vu e L^{PtY. □ 

Now, given / e ^ let us consider the equation 

- Am = / in n, Vm-n = 0 in 5^2. (5.13) 

We say that m e is a weak solution of (5.13) if 

J Vrn{x)Vpi(x)dx = ((/,</?)) V (p e (hi). (5.14) 


Lemma 5.4. Assume that q> d and let a e M. Then, there exists a unique weak solution 

of (5.13) satisfying that / mda; = a. Moreover, there exists a constant c> 0, independent 

Jn 

of{cL,f), such that for any F solving (5.12) we have that 

IVtuIli <c\\F\\l<!. (5.15) 

Sketch of the proof: Noticing that (5.14) is invariant if a constant is added to m, it suffices 
to prove the result for a = 0. Since q > d we have that q' <2 and so, by the Lax-Milgram 
theorem, existence and uniqueness for (5.14) holds in Using interpolation re¬ 

sults due to Stampacchia (see [41] and [42]), estimate (5.15) holds if Dirichlet-boundary 
conditions were considered (see e.g. [37, Theorem 7.1]). This argument yields the desired 
local regularity for m, which can be extended up to the boundary (which we recall that 
it is assumed to be regular) using classical reflexion arguments. 

□ 


Finally let us recall the following result about elliptic equations with measure data. 


Theorem 5.5 ([29], Theorem 1.2). Let f e 911(11). Then the unique solution u e lTQ^’^(r2) 
of the problem 


- Au = f in Q, w = 0 on 511 


(5.16) 


has the following regularity properties: 

(i) Vu e Vl/i^;"’^(ll)^ for all e e (0,1). 

j. 

(ii) More generally, Vu e ’ (H)'^, for all e e (0,cr(r)) where 1 < r < 

d - r(d - 1). 


d 

d-1 


and a(r) := 


Remark 5.2. We remark that the result about the uniqueness of the (renormalized) solu¬ 
tion of the problem (5.16) can be found in [f3]. Moreover, since the regularity results in 
Theorem 5.5 are local, these remain true if we use homogeneous Neumann boundary con¬ 
ditions instead of Dirichlet ones. In this context the solution is unique up to an additive 
constant. 
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